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ABSTRACT 

Same steady-state solutions are obtained for thermal and for shear- 
flow turbulence by expanding two-point nonlinear correlation equations in 
power series in the space variables. The correlation equations, which 
are for inhomogeneous turbulence, are constructed from the Navier-Stokes 
and energy equations. To make the problem determinate, the weak- 
turbulence approximation is used. Steady-state solutions are possible 
because of the presence of nonlinear production terms in the correlation 
equations. Because only the low-order terms are retained in the power 
series used in the expansions, the solutions are accurate only for small 
values of the space variables, and specific boundary conditions cannot be 
applied. The forms of the solutions show that critical values of parame- 
ters (similar to Rayleigh or Reynolds numbers) exist below which the tur- 
bulent fluctuations are zero. The main conclusion of the study is that 
the Navier-Stokes and energy equations (averaged for turbulent flow) can 
yield solutions in which the energy fed into a turbulent field by bupy- 
ancy or shear forces is equal to that dissipated by viscosity. 

INTRODUCTION 

Most of the analyses based on statistical turbulence theory have 
been made for a decaying turbulence that is initially generated by exter- 
nal means, as by flow through a grid. Sustained turbulence, such as 
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that produced by shear or buoyancy forces, has generally been investiga- 

4-7 

ted by using a phenomenological approach, or by using simplified equa- 

D 

tions in place of the equations for the real fluid. Statistical methods 
have also been used for sustained turbulence, 4 *"^ but the work is gener- 
ally limited to one-point correlation equations that by themselves do not 
lead to solutions, although they are of considerable schematic value. 

The work done to date offers little evidence that the Havier- Stokes 
and energy equations either are, or are not, capable of yielding solu- 
tions for steady- state turbulence. Studies of the effect of a uniform 
velocity gradient 1 ® and of a uniform temperature gradient and body 
force 11 on an initially isotropic turbulent field indicate that in those 
cases, although energy can be fed into a turbulent field by shear or 
buoyancy forces, the energy fed in is less than that dissipated, and the 
turbulence decays with time. If it were not for the abundance of sus- 
tained turbulent flows in nature, one might, on the basis of the available 
solutions, be led to the conclusion that steady-state turbulence will not 
occur. 

The present study is an attempt to provide evidence that the Uavier- 
Stokes and energy equations in averaged form can yield solutions for 
steady-state turbulence. The work is based on generalized two-point cor- 
relation equations that are constructed from the Navier-Stokes and energy 
equations by methods similar to those used by von K&rmdn and Howarth for 
isotropic turbulence.^ To make the problem determinate, the weak turbu- 
lence approximation (triple correlations neglected) is used. This approx- 
imation was also used by von Karman and Howarth for the case of low 
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Reynolds number turbulence. Although the approximation might be con- 
sidered somewhat restrictive, it appears to be the only reasonable basis 

of analysis, unless three- or four-point correlation equations are con- 
3 

sldered. Moreover, since we are studying sustained turbulence, we are 
more interested in the production terms in the equations than in the 
transfer terms, (in the case of shear-flow turbulence the velocity gra- 
dient causes energy transfer between wave numbers even when triple cor- 
relations are neglected.^®) 

12 

Sustained turbulence is essentially a nonlinear phenomenon. The 
nonlinear character of the two-point correlation equations is made evi- 
dent when the mean temperatures or velocities are eliminated by intro- 
ducing one point correlation equations into the twcwpoint correlation 
equations. Plane heat transfer and shear layers are considered. The 
correlation equations are expanded in power series in the space varia- 
bles to obtain algebriac expressions for the correlations. Because only 
the low- order terms are retained in the series, the solutions are accu- 
rate only for small, values of the space variables. 

The case of sustained thermal turbulence will be considered in the 
next section, after which sustained shear-flow turbulence will be taken 
up. 

SUSTAINED THERMAL TURBULENCE 

The term thermal turbulence as used here designates turbulence that 
is sustained by buoyance forces arising from temperature gradients and a 
body force* For steady-state turbulence, the energy fed into the turbulent 
field by buoyancy forces just balances that dissipated by viscous action. 



Correlation Equations 


Two-point correlation equations for homogeneous turbulence with a 
body force and a temperature gradient are constructed from the Eavier- 
Stokes aryj energy equations in reference 11* Here only the modifications 
necessary for inhomogeneous turbulence will be considered* From refer- 
ence 11 the correlation between velocity components at two points is 
given by 

dt J dx k 3x k P 3x^ J P dxj dx k dx k 
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■where the overbars Indicate averaged values* The subscripts can take on 
the values 1 , 2, or 3, »nd a repeated subscript in a term indicates a 
summation* The quantity Uj_ is an instantaneous velocity component, x^ 
and x^ are space coordinates at the points P and P 1 , t is the time, 
P is the density, v is the kinematic viscosity, p is the instantaneous 
pressure, is a component of Hie body force, *r is the fluctuating 
part of the instantaneous temperature, and 3 is the thermal expansion 
coefficient defined by 3 = - (l/p)(Sp/ST)p* For inhomogeneous turbu- 
lence it is convenient to introduce the variables r k s x k - x k a nd 
(x k ) m ^ 1/2 (x k + xi) ( see Fig. l)« Then Eq. (l) becomes 
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where the following transformations were used: 
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Equation (2) reduces to the ESnr^n-Hcwarth equation 1 if the turbulence is 

homogeneous and body forces are absent. 

In a similar way the following equations are obtained. 

Pres sure- velocity correlations: 


S _ 1 5 + JL 

dx£ 2 S(x k ) m Sr k 

S 2 




d 2 pu' 


4 ^C^kW^k ^ r k^ r kJ 


, i , i a Sv.i 

2 S(x z )u^r k 2 d(x k )n£r 2 dr 2 Sr k 



1 

1 


M 

SO 

^ ^(x 2 ) m S (x k ) m 

'Svj 

- pgk 

^ Stuj 

Sr 2 Sr k 

2 d(x k ) m dr k 


( 3 ) 


Temperature-velocity correlations: 
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Temperature-pressure correlations: 
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where T is the mean temperature and. a is the thermal diffusivity. 
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These equations reduce to Eqs* (7), (9), (15), and (19) of reference 11 
if the turbulence is homogeneous. It is also of interest to compare 
them with Eqs. (5), (7), and (8) of reference 10. 

Assume now that the only nonzero component of "g is in the negative 
vertical direction, and let 


g - - S3 


(7) 
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Also, let the temperature gradient be in the vertical direction so that 
it is given by ST/Sx^. The turbulence can then be homogeneous in hori- 
zontal planes, and vertical axes will be axes of symmetry. Let 

rf + rf - (8) 

If the turbulence is weak enough for triple correlations to be neglected 
and if we let i = j = 3 in Eqs. (2) to (6), 
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In order to eliminate the mean temperature gradient from the pre- 
ceding two-point correlation equations and to emphasize their no nlin ear 


character, we obtain a one-point correlation equation by sbustituting 
Eq. (5) into Eq. (ll) of reference 11 and by averaging. This gives 

ST . S 2 T 

+ — a 
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In the remainder of the section we will be concerned only with the steady- 

state case. Also, since the correlations change only in the x di- 

o 

rection, Eq. (14) becomes 
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where q 3 is the heat transfer per unit area and is independent of posi- 
tion. The quantity k Is the thermal conductivity. The temperature 
gradient at point P is then given by 


ST _ q 3 

Sx 3 k 

and that at point P 1 is given by 

_ ST’ = ^3 . 
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Substitution of Eqs. (16) and (17) into the two-point correlation 
Eqs. (9), (10), and (ll) gives, for the steady-state case, 
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Eqs. (18, (19), .and (20) together with Eqs. ('12 ) and (13) form a deter- 
minate set. (Note that *ru 3 and t'u 3 are special cases of ru 3 . ) 


Inspection of the equations shows that one possible solution is given 

when the correlations are all zero. In that case no turbulence will 

b 

occur and the heat transfer will be entirely by conduction. We would 
expect that if the fluid is heated from below (positive q 3 ) there will, 
in addition, by a nonzero solution, inasmuch as experiment indicates 
that turbulence can be set up for that case. Also, the presence of the 
nonlinear terms u 3 uj tu 3 , TU 3 (-r) tu 3 , and tu 3 t’u 3 in Eqs. (19) 
and (20) leads us to suspect that nonzero solutions exist. If those 
terms were not present, the no-turbulence solution would in general, be 
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the only pertinent solution; the equations in that case would he linear 
and homogeneous. In the next section the possibility of turbulent solu- 
tions of the steady- state correlation equations will be investigated by 
expanding them in power series. Before doing that, however , it is con- 
venient to convert them to dimensionless form by introducing the follow- 
ing dimensionless variables: 
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where A is a typical microscale of the turbulence. The microscale is 
used as a length because it is defined in terms of the shape of the cor- 
relation curve for small values of the space variables, and the solutions 
to be obtained are accurate only for small values of those variables. 

The quantity N t is a determining parameter for the thermal turbulence 
and is somewhat similar to a Rayleigh number. Eqs. (18) to (20) and 
Eqs. (12) and (13) become, in dimensionless form. 
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Series Expansions 

The plan now is to expand each correlation in Eqs. (2l) to (25) in 
a power series in |*, to substitute the series into the partial differen- 
tial equations, and to equate the coefficient of each power of |* to 
zero. This gives a set of partial differential equations that does not 
contain | . Each dependent variable in those equations is then ex- 
panded in a series in r^ to obtain ordinary differential equations that 
do not contain r^. Finally, expansion in (x 3 )* eliminates that varia- 
ble, and we end up with a set of algebraic equations that can be solved 
simultaneously to obtain values for the correlations. 
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In the present analysis we will consider only equations obtained by- 

setting the coefficients of g*®, rt^, r*"S (x^)*^, and (x^)*^ equal to 

zero. Thus the solution obtained will be accurate only for small values 

of those variables and will approach an exact solution only as £ , rj, ' 

and (x^)* approach zero. We will not be able to accurately apply 

boundary conditions that state, for instance, that a correlation is zero 

at given points in £*, r^, (x^)* space (e.g. , at walls or at «), 

since, in general, g*, etc. will not be small at the points the boundary 

conditions are applied. In lieu of boundary conditions we will introduce 

microscales that depend only on the shapes of the correlation curves near 

their origins. As defined by Taylor, a microscale is the distance at 

which the inscribed parabola at the origin of a correlation curve goes to 

zero. 4 For the law Reynolds numbers considered here (final period for 

decaying turbulence), the microscale of the turbulence differs but 

3 

slightly from the macroscale. The microscales used here are, in some 
cases, slight generalizations of the usual concept inasmuch as we include 
microscales associated with (xj^ as well as with £ and rj. We will 
also consider microscales for the case where the slope of the correlation 
curve at its origin is not zero. In the latter case a third degree rather 
than a second degree curve is inscribed in the correlation curve at its 
origin. 

The various correlations will, in general, have different microscales. 
For the sake of definiteness, they will here be arbitrarily taken as equal 
in most cases and will be designated by A. (The A used here is twice 
the usual microscale. ) In order to obtain the actual relation between 
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the microscales for specific boundary conditions (-which we are not spec- 
ifying here), it would he necessary to consider higher-order expansions 
of the correlation equations. The present analysis, however, should he 
adequate for determining whether or not reasonable solutions of the cor- 
relation equations exist for steady- state turbulence. Let 
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where the barred quantities in parenthesis are independent of £ . Then 
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(ujU^Jq, for instance, is the value of (U3U3 ) for £* = 0, and (^UjJg 
is (l/2)d 2 u 3 u^/d !* 2 evaluated at £* « 0. The odd powers of £* are 
omitted in these expressions because of symmetry. By using a Taylor 


series we can write TU 3 and t U 3 in Eqs. ( 22 ) and (23) as (Fig. l) 
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where (tu^) qo is the value of ru^ at (x 3 ) for £* = r 3 = 0. Sub- 
stituting Eqs. (26) to (28) in Eqs. ( 2 l) to (25) and setting the coef- 
*0 

ficient of £ equal to zero in each equation give 
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As mentioned previously, in order to obtain definite results we arbi- 
trarily assume that the microscales for the correlations are equal to A. 
Thus, in Eqs. (26) let u 3 u 3 = pu 3 = xu 3 = tt' = xp* = 0 for 
I* = l/A = 1. (Note that the actual correlations are not zero at those 
points, but only the inscribed parabolas given by Eqs. (26) i) Then 
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Substituting Eqs. (34) into Eqs. (29) to (33) and setting the coefficient 
of rtP and of in each equation equal to zero give (with 

(U 3 U 3 ) = -(ugu^) , etc. ) a set of ordinary differential equations in 

(xj)*. In Eqs. (34) we let ( u^u^ )^ , (puj)^, etc. equal 0 for r^ m ±1 
to obtain (^J ) Q2 - (P«p 
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° m 


(x-)*® equal to zero give (with (u*u 1) = -(u~u 1) > etc.). 


m 


3 3' 


02 


3 3' 


00 



- 17 - 


0 = -(pui) + 2(pu') + (ugii* ) - 12(u,ui ) + 2(mi) (36) 

5 001 5 010 002 3 3 000 3 000 


° = (u 3 u£) 


000 Pr 


^-Pr^ooo “ 2 ( ^ ) 00l‘ (Tp,) dl0 + 2 + Pr ) 


+ ( 1 -ff) ( ^on- 6 fr + 1 ) ( ^?ooo +(WT) oa 

°' 2( ^>ooo [*? ‘ * ( ^ooo] [<~W-^(~)ooo] 


| (pul) 

6 3 002 


(pul) - 6 ( pul ) 
3 Oil 3 ( 


-(tu’) 

3 010 


(tP T ) 0 02 + (^)oil - 6 ^, 


0= “I (tp*) 0 n + 2 ^c 


7T (pul) + 
2 3 012 




2(PU 3) 001 - “(pojloio -5 (™3> 011 + i) 


2 (^012 ~ 2 ^ t P*)o 01 “ 1C) ( T P ' ) oio 
( ^3^10 = (TT,) 010 “ C 


-2( TT 1 )q00 


Setting the coefficients of (x*)* 1 equal to zero gives 

° m 

f <^>001 - 2 ( S J ; 012 - 6( ^>001 - '^002 - <^>011 


3 002 


2 + 2 ( t ^')o12 “ 6(t£>' ) ooi = ( TT ' )o02 

2 ( pu 3^013 + 4 ( pu 3^002 " • lo (p u 3)o 11 = ( T1 *3)oi2 

I ( ^013 * 4( ^002 - 10 ^) 0 11 " 0 


- 18 - 


0 - PrfujuJl^fTuii^ - ( TP ’ ) 01E + 2(rp 1 ) 001 + § (l + Pf)< tu 3> 0 13 


4 ( 4 - * 10 fe + 1 ) (1 ^ ) 


Oil 


®on = (TT,) on = 0 


(49) 

(50) 


To set the microscales associated with (x3) m equal to A, let (u^u^)^, 

(p5J) , etc. in Eqs. (35) equal zero when (*3)* = ±1/2. (We use 

(x,)* = ±l/2 instead of ±1 as in the other cases, because of the differ- 

ence between the definition of (x*) and that of r 3 or £. See 
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B ' v 3 3' 0 02 v 3 3' 000 ' v -^3'o02 3 000 

(pu^) * -4(puj)^^, etc. With these relations, Eqs. (36) to (50) form 

a determinate set of algebraic equations that can be solved simultaneously. 
Equations (40), (48), and (41) show that 


to* >000 ” to ' )qii ~ toi). 


010 


(51) 


From Eqs. (42), (45), and (36), 

<^>ooo - h (™J>ooo • ife 

Combining Eqs. (43), (46), and (49) gives 


011 


(52) 




Pr 


3 Oil 4 _1_ 
Pr 


(u,uJL) (ru') 

o 3'qoO 3 


000 


+ A- (TT T )mn + 


(i-£) 




From Eq. (38) 


< TT '>000 ' 8 ^^ooo 


> 000 ' 1 

^ + 1 


W t 

■5— - Pt(tu’) 

Ft 3 000 


3' 000 


(53) 


(54) 
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Finally, Eqs. (43), (46), (52) to (54) and (37) combine to give 
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3 1 Pr - 1 \ 

32 " 128 Er + l/“ 


51^84 “ 2 + h ^ 


- 8 


Pr +1 v 
Ft 1 


i 1 pJ-L . Jl. ** - M 

4 Er + 1 ^52 128 Er + 1 / k 3' 000 . 


1 Pr 3 (Pr - 1 ) 
65,536 (Pp + -^2 



CD + 


000 



y *+ 


1 

128 


I 


Pr - 1) 
Pr + l) 


03)T 


where 


K+ 

= p~ - Pt(tu’) 

Er 3 y 000 


(55) 


and 

T “ 1 " 5l2 (T7 Pr) ^ TU 3^000 

Note that the original (ml) in the linear terms of Eq. (37) (that 

° 000 

equation having come from Eq. (10), the equation for tu^ ) canceled out 

of Eq. (55). However, Eq. (55) still yields values for (rui) , which 

* 000 

is the value of tuj for g " r 3 " ( x 5 ) ffl a because the original equa- 
tions were nonlinear in (tuJ) 

K 3 '000 

Positive values of correspond to negative temperature gradients 

(heating from below) , and turbulent solutions should exist for suffi- 
ciently high values of N^. Equation (55) was solved for Pr - 0.7 
(corresponding approximately to gases) and for several values of N^. 

As increases, (tu^), as well as the other correlations, becomes 

positive at Hj. = 88.1. Thus that value of Hj. can be thought of as a 
critical N-j- above which turbulence might occur. For = 100, 

(rui) j= 26.4, and for IL = 200, (ruJ) = 251. From Eqs. (52) to 
5 000 ^ 3 000 

(54), for = 100, (U 3 U 5 ) « 2.49 and (tt* )qoo = 287, and for 

vJvJvy 
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H|. = 200, ( u 3 u 3 )qqq 30 27.3 and ( t t’)ooo “ 2,414. Thus the correlations 
have the correct signs and the correct trends as Hj. increases above its 
critical value. For values of belov the critical value, including 

negative Sj. (positive temperature gradient), the correlations become 
negative according to Eqs. (52) to (54), so that those equations do not 
yield possible solutions. In that case we should use the other solution 
of Eqs. (21) to (25), that is, the no-turbulence solution. Note that 
above the critical both the turbulent and the nonturbulent solutions 

are possible according to the equations j consequently, the fluid is not 
necessarily turbulent for all values of N-j. above the critical value 
considered here. 

One other quantity that should be considered is the eddy diffusivity 
for heat transfer e h , -which is defined as 


TUj TUj 


(56) 


03 

k 


™3 

a 


In dimensionless form, at (x,)* = 0, 

° m 


<™ 3 > 


000 


'h N. 

|| - Et(tu^) 


000 


From the preceding computed results, is positive and increases as 
N^. increases above its critical value so that the values obtained for 
e h are also reasonable. 

It might again be mentioned that the preceding solutions are for the 
case where the microscales are all equal. Thus the results probably do 
not correspond numerically to those for particular boundary conditions. 
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They are given here to indicate that reasonable solutions for steady- 

state turbulence can he obtained from the correlation equations. 

Equation (55) is quadratic in (tuj) and even for unequal micro- 

scales the solution for (ruT) is of the form 

^000 

( ™?ooo * a + ™t * + eS t + f (57) 

where a, b, etc. are functions of Prandtl number and the microscale 
ratios. (The microscale ratios might in turn be functions of but 

those functions are probably slcrwly varying. ) Setting (tuJ) =0 in 
Eq. (57) gives two values for a critical (for given values of a, b , 
etc.), but, of course, only one of those would be expected to be physically 
realizable. 

SUSTAINED SHEAR-FLOW TURBULENCE 


Correlation Equations 

General two-point correlation equations for incompressible turbulent 
shear flow were obtained from the Navier-Stokes equations in reference 10 

as 


g dUi 

?t U i U J +U k U J 3^ + u i u k 




+ 1 K u K + ^i^j) + ( u i u j u k - u i u k u j) 


P 1 2 
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3GSE PU J + UiP J + ^I UiP ' P ^ 


+ 1 v d u i u .i 

+ 2 v 3Tx k ') S^y ■ 
^ m m 


a ujuj 


+ 2v 


(58) 
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1 Su i u k , du i u k 
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4 d ( x k) m 5 x k) m + ^k^k + ^k^k 
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1 uu i u l u k 
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( 59 ) 
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9t 
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1 5 u z u k uj t 1 5 u Z u k u j , 1 ^JWjS ^ U * U k U i 

" 4 2 2 d ( x k) m ^Z " 


(60) 


where Uj_ is the fluctuating part of an instantaneous velocity component, 
Ui is a mean velocity component, and the other quantities have the same 
meanings as in the equations in the preceding sections. As before, the 
unprimed quantities are measured at point P and the primed ones at P 1 . 
The vector configuration in Fig. 1 applies to the present equations. A 
one-point equation given in reference 10 is 


dU. 


^ - - p & + sr r ss: ‘ ^ 


(61) 


It is suggested in reference 10 that Eqs. (58) to (61), together with 

* V 

higher -order equations , should constitute a solution to the turbulent 
shear-floir problem. Further evidence that this is the case is given 
herein. 

In this analysis it is assumed that the mean velocity is in the x^ 


direction and that changes in the mean quantities can take place only in 
the Xg direction; consequently, *we have plane shear layers. If the 


turbulence 1 b ve&k enough for triple correlations to be neglected and if 
-we let i = j = 2 in Eg.. (58), that equation becomes 

i/il 


St U 2 U £ " ’( U l " U l) ^ *2*2 - ^ £(x 2 ) m + 


. d ,r-rw ^ 3 r-rl , 1 B U 2 U 2 . 0 u U 2 U : 

+ (p^K-r) - Pu 2 > + 2 v S( X2 ) m 5;x 2 ) m + 2v ^ 


S U 2 U^ 
r k 


(62) 


Similarly, if i ■ 1 and j « 2 in Eqs. (58) to (60), -we obtain 

dU-, 


a r r 1 TT x a r 1 

St U 1 U 2 “ “ U 2 U 2 ^ “ ( U 1 - U l) ^ U 1 U 2 ~ p 


1 a — , 

2 d ( x 2) m UlP 
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(63) 
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au x au 2 u 2 
axg ar^ 


(64) 


and 


1 

P 


a Uj^p' 


a u x p’ 


w d u x p 

* + + ^ 


-2 


au^ au x u^ 


(65) 


The one-point equation, Eq. (61), can be used to eliminate the mean ve- 
locities from these equations. If the turbulence is steady state, 

Eq. (61) becomes (for the plane shear layers considered here vith no mean 
pressure gradient) 
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or 


dU- 


v 357 ' u i u 2 = “ 


_T 

P 


( 66 ) 


where t is the total shear stress (made up of the laminar and the tur- 
bulent shear stress) and is independent of position. Then the velocity 
gradient at point P is 

3u n 


= jr_ . 1 
c 2 “ pv 


\ = — — + — U-i Up 

OXp pv V X 6 


(67) 


and that at P* is 


At a general point P , 


SD i jr_ . 1 — r-r 

dxg “ pv v 


T 


( 68 ) 


(69) 


Integrating Eq. (69) between P and P 1 , that is, from (x 2 ) m - 1 * 2/2 to 
(*2) m + T z/ 2 ( Fi S- !)> gives 

’< x 2) m +r 2/ 2 


Un - U-, = — ro + - 

X ± py C V 


u l u 2 ^2 


(70) 


( x 2) m - r 2/ 2 

The substitution of Eqs. (67), (68), and- (70) into Eqs. (62) to (65) 
gives, for steady-state turbulence, 
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£7 r 2 



u l"2 ^2 ] ^ 2 U 2 
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d Zu 2 u 2 
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(71) 
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* 



and 


1 

P 


<z r 

a u n p' 


-s2 — 

a u. 


4 a(x 2 )| + ^(x 2 ) m ^ 2 + 


■' 2 (^ + 7 u i u 2)Jj“l u 2 ( 74 > 


Equations (71) to (74) form a determinate set. A possible solution 
for these equations is that all the correlations are zero, as "was the 
case for the equations for thermal turbulence. In that case no turbulence 
exists and the shear stress is produced by molecular action. We are here 
mainly interested in turbulent solutions that might be possible because 
of the nonlinear terms in the equations. Of these nonlinear terms the 
first term of Eq. (71) and the second term of Eq. (72) have been inter- 
preted' 1 '® as Fourier transforms of transfer terms that transfer energy 
between eddies of various sizes. The remainder of the nonlinear terms, 
that is, the first term on the right side of Eq. (72) and the last terms 
in Eqs. (73) and (74), might be interpreted as production terms. Some of 
the nonlinear terms here are more complicated than those for thermal tur- 
bulence, where they were all simple products of correlations. The equa- 
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tions can be -written in dimensionless form by introducing the following 
dimensionless -variables: 


* A UjuI 

u i u j --r 1 ' 


H . IL- r * = — fx ) 

» pv 2> r i A ' ^ X 2 J 


(x 2 } 

* m 
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A * pu 


'J 


pv 3 


where the microscale A is again used as a length scale because the 
solutions to be obtained are accurate o nly for small values of the space 
variables. The parameter N s is a determining parameter for the shear- 
flow turbulence and has some similarity to the square of a Reynolds num- 
ber. Equations (71) to (74) became, in dimensionless form, 
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(76) 
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Series Expansions 

The expansion in pover series of the correlation equations for shear 
flow is similar to that for thermal turbulence. Most of the discussion 
at the beginning of the section "Series Expansions" for sustained thermal 
turbulence therefore applies also to shear-flow turbulence. The main 
difference between the two cases is that we had the three independent 
variables g, r^, and (xj)^ for thermal turbulence, whereas for shear 
flow we have the four variables rp r 2 , rj, and (x£) m - This is because 
there is no axial symmetry for shear-flow turbulence as there was for 
thermal turbulence. The final expanded correlations will then have four 
subscripts instead of three. As before, all of the microscales will at 
the beginning, for the sake of definiteness, be assumed equal to A. In 
the present case, however, it will be found necessary to later modify 
that assumption to obtain reasonable results. 

We will first obtain expressions for the nonlinear tenns in Eqs. 

(75) to (78) in terms of the expanded correlations. The quantities 
u l u 2* 8111(1 u-[u£ X can be expanded in a Taylor series to give, for 
small rS, 

c r* 

U 1 U 2 = (^>000 " ~2 d ( x 2 )* ( U ] - U 2 ) 000 


( 79 ) 
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and 


^2 " ^l^QOO + 2 d(x 2 ) m W U 2^ 


OOO 


(80) 


where (uqiip 000 is the value of u-jUg at the point (xg)* with 

r 3 “ r 2 * r l “ ®* 1116 quantity (u-jU^) , in turn, can he expanded as 

'V^ooo - <V?oooo + ( V5W**>: 2 < 81 > 

where (u^ujj^^ is the value of for (x g )* = r* = r* = r* = 0. 

The odd powers of (x ? )* are omitted to make (uTuT) symmetric about 

c m - 1 - ^ 000 


(x 2 )* = 0. To relate ( u i u 2)oo 02 to ( u i u 2)oooo' set the microscale for 

(vDoooo 


0 for 


and ( x 2 ) m equal to A by letting ( u i u 2)qoo 
( x 2^m = V 2 - This gives (u^y$) 0002 - -^u^Jqqqq. Then Eqs. (79) 
to (81) yield 
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Let 
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u 2 u£ = (u 2 u*) 0 + (u 2 u^) 2 r* 2 

- ( ”?Po + ( V? } 2 T f 
- ( i ^) 0 + ( i 5 p 2 r f 

U 3 P 7 * * (^aJ r ) 0 + ( u iP * ) 2 r I' 


where the quantities in parenthesis are independent of rtj, and (uguj,)^, 
for instance, is the value of UgU^ for r^ = 0. The odd powers of r^ 
are omitted because of symmetry. Substituting Eqs. (82), (83), (85), and 
(86) into Eqs. (75) to (78) and equating to zero the coefficient of rg 
in each equation give 

° - - |n b + ( : ^2 ) 000 o]r| ^5: (^J) 0 - | S (x 2 )* 
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In order to eliminate rj> from the correlation equations, let 

( u 2 u <Po = ^ U 2 U 2^oo + r 2 + ^ U 2 U 2^02 + ^ U 2 U 2^ 0 3 r % 
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(u l p,) o " (u l p,) 00 + (U l p,) 01 *2 + (u l p ’) 0 2 r 2 2 + (u l p,) 03 *2 


*3 


J 


where, for instance, ( UgUg ) qq is the value of ujju*? ^ or r 3 = r 2 = ®* 
Substituting Eqs . (91) into (87) to (90) and equating v to zero the coeffi- 
cient of t?P give* with (ugUg)^ * -(ugUg)^, etc. , a set of partial dif- 
ferential equations with independent variables r* and (x^)*. Setting 


the microscales for ( UgUg ) q , ( u i u 2 ) q ’ ( P u 2 ) q 3 an< * ( U ]_P ' ) q in Eqs. (91) 


equal to A by letting (ugUg)^, (nin^o* e ^ c " “ ® ^ or r 2 = we 
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To eliminate r* 

^oo-^ooo* 

(^ ) 0l“ ( ^ I 2 ) 010 + 
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from the correlation equations, let 
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(92) 


(-1P')01 = ^010 + ^011 r l + ( u lP’) 0 12 r l 2 + (^lP r )oi3 r l 3 

Substituting Eqs. (92) into the correlation equations with independent 
variables rj and (xg)* (not shown) and equating to zero the coeffi- 
cients of r^° and r^ 1 give (withfugU^)^ * -(ugUg)^, efcc ') a set 

ordinary differential equations in (x ? )*. 

" m 

Finally, to eliminate (xg)*, let 
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justing the coefficientB of (xg)* P to zero gives 
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Equations (106) and (104) give 


(U 2 U 2 J 0011 * (U 2 U 2*0101 * ° 


( 110 ) 


since (ugii^)^^ and (u 2 u£) 01Q3 are proportional respectively to 


( U 2 U 2 ) ooil and (u 2 u 2) 010 i' 


Equation (94) emphasizes the importance of the pressure-velocity 
correlations for shear-flov turbulence. In the absence of those correla- 
tions (^^^qqqq 'would be zero (since (ugUg)^^ is proportional to 

-(^Poooo) and ttiere "would be no turbulence. This is not surprising, 
since the turbulent energy is fed into the turbulent field through the 
u^u^ component and is distributed between the various components of the 
energy by the pressure-velocity correlation terms. 
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As in the case of thermal turbulence , it is the presence of the non- 
linear or quadratic terms in some of Eqs. (94) to (109) that makes possible 
a nonzero solution of those equations. In the absence of those terms the 
equations would be linear and homogeneous and, in general, would have only 
a no-turbulence solution. 

To relate (u^up 00Q2 to (ugu^)^^, etc. let the microscales asso- 
ciated with (x 9 )* in Eqs. (93) be equal to A by letting, far instance, 

^ m 

'V'Pooo * 0 for - ±1/2 - 11118 glT “ 'V'Pocos * - 4( V?oooo’ 

(^2>0003 ‘ - 4 <^0001> ''^F'oOlZ " '^“^OOIO’ etC ' Dsln ® these 
relations we get, when substituting Eqs. (96), (98), (100), (102), (108), 

and (110) into (94), 
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Note that ( U 2 U 2)()000 cance ^- s ou ‘ t this equation leaving a quadratic 

equation in (u-ju|) 00 . Thus we get the somewhat unexpected result that 

Eq. (94), which was originally obtained from Eq. (62) (the equation for 

UgUg ) > gives, when combined with the other equations, a solution for 

Similarly, Eq. (95) will give a 

Solution of 
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solution for (u 0 u’) rather than for (u_ui) 
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Eq. (ill) gives 


(u,ui) = -2H g ± 

1 * 0000 B 
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( 112 ) 


This equation can give a negative value of u^u^ for positive E s and 

a positive u u ’ for negative H g , as it should . However, the critical 

No (the value of KL for (uTuT) = 0) is imaginary. Moreover, if we 

B s 1 2 0000 

calculate the eddy diffusivity for momentum transfer from 
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which in dimensionless form at (xj* becomes 

c m 
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0000 


N b + (u^uj) 


(113) 


we find that Eq. (112) will not give a positive solution for e as re- 
quired. This only indicates that all the microscales cannot be equal as 
assumed. If, for instance, we take the microscale for (pul) in 

Eqs. (93) as A/2 instead of A, we get (p'ug)^^ = -16(pu£) 010Q and 
(pug) 0103 - -16(pu^) o ‘ the o'ther microscales are taken equal to A, 


<vP, 


- -2N S ± a/E“ + 17,024 


where the positive sign is taken for positive N g and the negative one 

for negative N s . Equation (114) gives a critical value of N s (at which 

(u-^Ug) = 0) of ±75.3. Also, substitution of Eq. (114) into (113) gives 

a positive eddy diffusivity, which has the correct trends with increasing 

N s . Physically realizable results for ( u-, uA ) and € were also 

± ^0000 

obtained by letting the microscales for ( pug ) ^ , (pug)^, and (pug)^^ in 
Eqs. (91), (92), and (93) be equal to 4A, the other microscales being set 
equal to A as before. 

To obtain (ugUg) as a function of N s , we solve Eqs. (103) to 

(llO) and Eq. (114) simultaneously. The values for microscales used for 
obtaining Eq. (114) were used here. For (ugUg)^^ “ N s * ±75.3, as 
was obtained for ( u i u 2^o000 “ 0 in (U4). For N g = ±100, 

( Vp0000 " 93 ' *** for N s - ±150 > (^Poooo “ 197 - ^ ( VpOOOO' 
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as veil as and e, has the correct trends and signs. Note 

that possible solutions are obtained for either positive or negative N s , 
whereas for the case of thermal turbulence solutions vere obtained only 
for positive Nf This is to be expected, since the turbulence should not 
be affected by the sign of the shear stress. For thermal turbulence, on 
the other hand, the turbulence should exist only for heat transfer in the 
positive direction (negative temperature gradients). 

As in the case of thermal turbulence, the results obtained in the 
preceding paragraphs probably do not correspond numerically to those for 
particular boundary conditions because the microscale ratios vould be 
different. They are given to shov that reasonable results can be ob- 
tained from the correlation equations. In the general case, Eq. (ill) 
can be -written as 
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vhere A and a are functions of the microscale ratios (-which might in 

ve 


turn be weak functions of N e ). Solving Eq. (115) for (u-^Ug) 
get 
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For 


0000 

(u^Ug)gQQQ = 0, N B = Thus A and a should both be either 

positive or negative, at least in the vicinity of ( u-^Ug ) qqqq “ 0. Com- 
parison of Eq. (116) to (57) shovs that the expressions for both shear- 
flov and thermal turbulence are solutions of quadratic equations, although 
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the forms differ somewhat. 


CONCLUDING REMARKS 


By expanding the two-point nonlinear correlation equations in power 
series, reasonable solutions were obtained for Doth thermal and shear- 
flow turbulence. By reasonable it is meant that the correlations and 
eddy diffusivities had the correct signs and trends. Moreover, critical 
values of the determining parameters, below which unphysical turbulent 
solutions occurred, were obtained. Below the critical values the no- 
turbulence solution of the correlation equations was therefore appro- 
priate. Above the critical values the equations showed that the fluid 
could be either turbulent or nonturbulent. For shear-flow turbulence the 
same solutions are obtained for either positive or negative shear stress, 
whereas for thermal turbulence it was necessary for the heat transfer to 
be positive (negative temperature gradient). These results are, of course, 
to be expected if the equations yield reasonable solutions. The results 
obtained would not be expected to correspond numerically to those for 
particular boundary conditions inasmuch as the microscale ratios corre- 
sponding to those conditions were not determined. To determine those it 
would be necessary to carry higher order terms in the expansions and apply 
the particular boundary conditions of interest. The presence of pressure- 
velocity correlations in the equations was found to be indispensible if 
steady-state shear-flow turbulence is to exist. The pressure-velocity 
correlations were of less importance for thermal turbulence. 
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